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Abstract 

Heat  transfer  and  fluid  flow  through  porous  media  was  investigated  using  numerical  simulations 
and  experiment.  For  the  numerical  simulations,  two  models  were  created.  The  first  consisted  of  a 
two-dimensional  numerical  model  created  in  MathCAD  and  was  solved  using  the  finite  difference 
approach.  The  MathCAD  model’s  flow  in  the  porous  media  was  described  by  the  Brinkman- 
Forchheimer-extended  Darcy  equation.  The  second  model  consisted  of  a  computational  fluid 
dynamics  (CFD)  porous  media  model  using  Fluent™  and  was  solved  using  the  finite  volume 
approach.  Both  models  assumed  constant  fluid  phase  and  properties.  Pore  diameters  were  held 
constant  for  each  simulation;  two  different  porosities  were  investigated.  Boundary  conditions  were 
applied  at  the  wall  in  which  the  temperatures  of  the  fluid  and  the  porous  media  were  determined  by 
coupled  energy  equations.  The  effects  of  the  boundary  condition,  the  Reynolds  number,  porosity, 
and  heat  input  were  examined. 

The  experimental  investigation  consisted  of  a  flow  channel  with  a  porous  media  section  that  was 
heated  from  below  by  a  heat  source.  The  variation  of  temperature  of  the  fluid  in  the  porous  media 
was  measured  along  the  centerline  and  along  the  top  wall  and  bottom  wall.  The  heat  source 
temperature  and  the  fluid’s  inlet  and  outlet  temperatures  were  also  measured.  The  results  of  the 
numerical  and  CFD  models  as  compared  to  the  experimental  data  for  fluid  flow  through  porous 
media  are  presented  in  the  paper. 
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v  volumetric 

x  x-direction 

y  y-direction 

1.  Introduction 

When  fluid  flows  through  a  porous  media  such  as  a  heat  exchanger,  packed  bed  regenerator,  etc.  the 
solution  to  the  equations  for  fluid  flow  become  difficult  due  to  the  tortuous  geometry  in  which  the 
fluid  must  navigate.  In  addition,  the  changes  in  porosity  that  flow  can  encounter  add  to  the 
difficulty  in  obtaining  a  solution  for  the  governing  equations.  Through  the  years  many  researchers 
have  studied  the  flow  characteristics  and  have  added  terms  and  laws  to  the  equations  in  order  to 
facilitate  the  conditions  that  the  fluid  sees  as  it  flows  through  the  porous  media.  Darcy’s  Law,  the 
Brinkman  extension  and  Forchheimer  are  typical  additions.  Also,  when  dealing  with  the  heat 
transfer  in  the  porous  media,  many  researchers  have  derived  boundary  conditions  for  the  model, 
such  as  local  thennal  equilibrium  and  local  thennal  non-equilibrium.  Each  of  the  two  cases 
mentioned  above  require  their  own  boundary  conditions  to  be  satisfied  resulting  in  different 
solutions  for  the  governing  equations. 

Over  the  last  decade  there  has  been  an  increased  amount  of  research  into  the  use  of  porous  media  as 
a  heat  transfer  enhancement.  Porous  media  heat  exchangers  (PMHE)  offer  a  large  amount  heat 
transfer  to  occur  within  a  small  volume.  Due  to  the  large  amount  of  surface  area  within  the  porous 
media,  a  PMHE  can  have  a  large  volume  to  surface  area  ratio.  As  opposed  to  the  local  thennal 
equilibrium  model  (LTE),  which  assumes  that  the  solid-phase  temperature  is  the  same  at  any  point 
in  the  porous  media,  the  local  thennal  non-equilibrium  model  (LTNE)  assumes  a  finite  difference 
between  the  fluid  and  the  porous  media  and  solves  for  their  respective  temperatures  separately. 

Numerous  authors  have  studied  the  LTNE  model’s  application  to  PMHE.  Pavel  et  al.  [2004] 
examined  the  influence  different  diameter  porous  media  had  on  the  Nusselt  number.  Their 
numerical  code  lacked  accuracy  due  to  the  restrictions  of  the  boundary  conditions  not  applying  to 
the  physical  model.  Their  conclusion  was  that  porous  media  0.8  times  the  inner  diameter  of  the 
pipe  offered  the  best  enhancement  of  heat  transfer  while  not  significantly  increasing  the  pressure 
drop.  Alazmi  and  Vafai  [2002],  study  how  different  boundary  conditions  can  affect  the  heat 
transfer  model.  Different  parameters  such  as  the  porosity,  Darcy  number,  Reynolds  number,  inertia 
parameter,  particle  diameter,  and  solid-to-fluid  conductivity  ratio  were  analyzed.  Jiang’s  et  al. 
[2001]  provided  an  in  depth  numerical  model  of  heat  transfer  in  a  channel  packed  with  various 
material.  The  numerical  model’s  output  is  compared  with  experimental  data  for  verification.  The 
authors  test  various  boundary  conditions  and  settle  on  the  ideal  constant  wall  heat  flux  boundary 
condition,  which  describes  the  heat  flux  transferred  by  the  fluid  and  by  the  solid  phases  are  equal. 
Also  discussed  are  the  effects  of  viscous  dissipation  on  the  heat  transfer  coefficient,  which  was  done 
using  a  numerical  model.  These  and  further  considerations  of  thermal  non-equilibrium  as  well  as 
the  thermal  equilibrium  models  and  other  boundary  conditions  are  outlined  in  Alazami  and 
Vafai[2001],  Jiang  et  al.  [1999],  Khashan  and  Al-Amiri[2005],  Ochoa-Tapia  and  Whitaker[1997], 
Jeng  et  al.  [2006],  and  Jen  and  Yan  [2005]. 

2.  Numerical,  CFD  and  Laboratory  Models 

Heat  transfer  and  flow  through  porous  media  was  modelled  with  MathCAD  and  Fluent™.  The 
effects  of  the  boundary  condition,  the  Reynolds  number,  heat  input,  and  porosity  were  all  examined. 
The  numerical  model  assumed  constant  fluid  phase,  constant  pore  diameter  per  porosity  experiment, 
viscous  dissipation,  no  conduction  in  the  x-direction  through  the  porous  media,  and  thermal  non¬ 
equilibrium  in  which  the  temperature  of  the  fluid  and  the  porous  media  are  at  different 
temperatures.  The  numerical  model  was  two-dimensional  and  solved  using  the  finite  difference 
approach.  The  momentum  equations  were  solved  first  to  obtain  the  velocity  profiles.  These 
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velocity  profiles  were  then  used  in  the  energy  equation  to  solve  for  the  heat  transfer  characteristics 
of  the  fluid  and  the  porous  media. 

The  Fluent™  model  was  based  on  flow  through  porous  media  and  solved  using  the  finite  volume 
method.  The  Fluent™  model  had  a  constant  heat  flux  applied  to  one  side  of  the  porous  media  with 
fluid  travelling  through  the  porous  media.  The  permeability  and  inertia  coefficient  were  defined 
along  with  other  fluid  characteristics  within  the  Fluent™  model. 

The  experimental  set-up,  similar  to  Boomsma  and  Poulikakos  [2001]  and  Jiang  et  al.  [2001], 
consisted  of  a  flow  channel  with  a  porous  media  section  that  was  heated  from  below.  The 
temperatures  of  the  fluid  in  the  porous  media  were  measured  down  the  centerline  and  along  each  of 
the  walls.  The  heater  block  temperature  was  also  monitored,  as  well  as  the  fluid’s  inlet  and  outlet 
temperatures.  The  pressure  drop  was  measured  across  the  porous  media  section  in  order  to 
determine  the  variables  used  in  Fluent™  and  the  numerical  model.  Finally,  the  results  of  the 
numerical  and  Fluent™  models  were  compared  to  the  experimental  data. 

3.  Numerical  Two-Dimensional  Flow  Model  for  Porous  Media 


The  steady-state  two-dimensional  governing  equations  for  single  phase  fluid  flow  in  an  isotropic, 
homogenous,  porous  media  based  on  the  Brinkman-Darcy-Forchheimer  model  with  constant 
properties  has  been  presented  by  numerous  authors  and  can  be  written  as; 

— f  ( u  •  V)u  =  -VP  ~  —  u  +  -j=  \Vmadu  +  —  V2u 
£2  K  4k  £ 

Equation  (1)  was  separated  into  x  and  y-direction  components  and  represented  non-dimensionally 
as; 

For  the  x-direction; 
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Taking  the  curl  and  using  the  definition  of  the  velocity  in  tenns  of  the  stream-function  yields; 
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Equation  (4)  represents  the  stream  function-velocity  model  for  flow  through  porous  media. 
Equation  (4)  was  solved  by  an  iterative  approach,  as  presented  by  Gupta  et  al.  [2005]. 


In  order  to  solve  equation  (4)  it  was  necessary  to  obtain  the  permeability  and  inertia  coefficients. 
The  Forcheimer  equation  was  used  to  obtain  these  values 
pFL 
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By  substituting  in  C, 
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. _ equation  ( 1 )  can  be  rewritten  as: 
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A p  =  CxuD  +  C2uD  (6) 

Then,  by  obtaining  the  pressure  drop  at  different  Darcian  velocities  across  the  porous  media  and 
plotting  the  results  the  coefficients;  Cl  and  C,  were  detennined  by  a  second  order  polynomial  curve 
fit.  From  knowing  the  coefficients  it  was  possible  to  obtain  the  permeability,  K  and  the  inertia 
coefficient,  F. 
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These  values  of  the  permeability  and  the  inertia  coefficient  were  used  in  the  numerical  model  as 
well  as  the  Fluent™  models. 

After  knowing  the  permeability  the  Darcy,  Da,  number  could  be  found  for  the  porous  media  flow; 
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4.  Numerical  2-D  Heat  Transfer  Model 


The  fluid-phase  and  solid-phase  energy  equations  for  the  thennal  non-equilibrium  model  with 
consideration  of  thennal  dispersion,  viscous  dissipation,  constant  porosity  and  constant  properties, 
and  neglecting  conduction  in  the  porous  media  in  the  x-direction  due  to  the  high  Peclet  number 
flows  being  investigated  (l  18  <  Pe  <  416) ,  were  as  follows;  for  the  fluid  phase; 
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Non-dimensionalizing  (7)  yielded; 


8X 


1 


umsXf 


Da  Re  Cp  ,qwH 
And  (8)  yielded; 


v  dY 


1 

a 

f  skf  +  Xd  ^ 

QOf " 

^Re  Pr 

dY 

l  A/  J 

dY 

ul  F XX 


+ 


/ 


CpfqwH  jDa 


UmagU  v  + 


Nu 


sf 


Re  Pr 


d2e 

dY2 


(1  s) 


dY2 


(9) 


(10) 


The  equations  were  descretized  and  then  solved  using  the  implicit  finite  difference  method  and 
using  the  velocities  found  by  the  two-dimensional  stream  function-velocity  model  equation  (4). 

The  boundary  condition  at  the  heated  wall  for  the  thermal  non-equilibrium  model  was  the  same  as 
mentioned  by  Jiang  et  al. [2001]. 
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The  Reynolds  number  could  be  represented  by  the  following  correlations  Jiang  et  al.  [2004],  Alazmi 
and  Vafai  [2001],  and  Boomsma  and  Poulikakos  [2001]. 
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The  equation  for  the  Reynolds  number  was  used  in  this  research,  unless  otherwise  noted,  was 
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The  local  heat  transfer  coefficient  was  calculated  as; 
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The  local  temperature  of  the  heat  transfer  surface  was  calculated  using  the  measured  temperatures 
of  the  wall. 
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Other  equations  for  the  Nusselt  number  were 
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Jiang  et  al.  [1999]  presented  a  modified  thermal  dispersion  equation; 
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The  equations  were  used  along  with  the  data  obtained  from  the  laboratory  experiments  to  model  the 
heat  transfer  in  porous  media.  These  results  were  then  compared  to  the  output  from  the  Fluent™ 
model  and  from  the  numerical  model. 


5.  Fluent™  Model 

The  porous  media  model  in  fluent  was  a  two-dimensional  model  consisting  of  a  flow  channel  with  a 
mid-section  representing  the  porous  media.  Gambit  was  used  to  create  the  geometry  and  the  mesh 
of  the  model.  The  mesh  was  finer  at  the  entrance  and  exit  of  the  porous  media.  The  grid  spacing 
for  the  porous  media  was  decreased  in  order  to  increase  the  accuracy  of  the  calculation  used  to 
determine  the  flow  characteristics  and  heat  transfer  within  the  porous  media. 
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Figure  1.  Temperature  Plane  Representation  in  the  Porous  Media  Model 


The  dimensions  of  the  Fluent™  model  matched  the  dimensions  of  the  physical  model  used  in  the 
laboratory  experiments.  Figure  1  shows  the  planes  used  to  record  temperature  data  within  the 
porous  media.  The  planes  are  grey  in  color.  They  are  in  both  the  x  and  y-direction  and  are  located 
at  the  entrance,  exit,  top,  bottom,  horizontally  through  the  middle  as  well  as  vertically  through  the 
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middle.  The  horizontal  planes  represent  the  location  of  the  thennocouples  used  in  the  laboratory 
experiment.  The  flow  in  Figure  1  was  from  left  to  right. 

Fluent™’s  porous  media  model  contained  settings  for  the  following  variables;  the  viscous 
resistance,  the  inertia  resistance,  and  the  porosity.  The  viscous  resistance  and  the  inertia  resistance 
values  were  obtained  from  the  experiment. 


The  boundary  condition  used  by  Fluent™  for  the  heat  flux  was; 
q  =  l,f(T,-Tr) 

Whereas,  the  fluid  side  heat  transfer  was  calculated  using 

=  -T, 


rdT^ 


dy 


(20) 


(21) 


For  post  processing  Fluent™  was  used  as  well.  The  post  processing  involved  the  output  files  of 
velocity  in  the  x  and  y-direction,  the  heat  input,  and  the  temperatures.  Each  of  the  laboratory 
experiments  were  replicated  in  Fluent™.  The  data  was  collected  and  graphed  for  comparison. 


6.  Experimental  Apparatus  and  Data  Collection 


The  physical  model  used  in  the  experiments  was  very  similar  to  the  model  used  by  Jiang  et  al. 
[2001],  The  model  assumes  the  porous  media  was  isotropic  and  homogeneous;  the  fluid  was  single 
phase,  constant  porosity,  and  constant  properties.  The  lower  side  of  the  porous  media  received  a 
constant  heat  flux,  while  the  upper  side  was  adiabatic.  The  test  apparatus  is  shown  in  Figure  2.  The 
fluid  channel  dimensions  were  .0381m  x  .0381m  x  ,0127m.  The  size  of  the  heated  test  section  was 
.038m  x  ,038m.  The  heater  block  in  the  lower  plate  of  the  channel  received  a  constant  heat  flux,  qw . 
The  flow  entered  the  channel  with  an  average  velocity,  uin ,  and  constant  temperature,  Tin . 
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Figure  2.  Experimental  Apparatus 


The  test  section  was  made  from  lexan,  and  consisted  of  an  upper  and  lower  section.  The  lower 
section  contained  the  flow  channel  and  section  for  the  porous  media.  The  porous  media  was 
obtained  from  ERG  Aerospace,  Inc.  constructed  of  aluminum  6101  with  constant  porosities.  Two 
different  pore  sizes  were  tested;  10PPI  (pores  per  inch)  and  40PPI.  The  upper  section  also 
contained  the  inlet  and  outlet  holes  for  the  fluid  as  well  as  the  thermocouple  insertion  points.  The 
lower  section  contained  the  heat  block  constructed  of  38mm  x  38mm  x  19mm  thick  aluminum 
which  held  two  120V/20W  heaters  6.35mm  diameter  and  25.4mm  length. 


Equally  spaced  thermocouples  were  inserted  mid-way  into  the  porous  media  along  the  horizontal 
center  line  of  the  test  section,  along  the  wall  at  the  top  corner,  and  along  the  bottom  comer  in  order 
to  monitor  the  fluid’s  temperature  at  each  of  the  locations.  The  heater  block’s  temperature  was  also 
monitored  with  a  thennocouple.  The  inlet  and  outlet  temperature  of  the  fluid  was  measured  with  a 
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thermocouple  placed  at  the  inlet  region  and  exit  region  of  the  lexan  test  apparatus.  Inlet  pressure 
was  monitored  using  a  pressure  transducer.  The  mass  flow  rate  of  the  water  was  controlled  by  3- 
piston  pump.  For  each  experiment,  water  was  used  as  the  working  fluid,  passing  it  through  the 
porous  media  while  holding  fixed  the  flow  rate,  inlet  temperature,  and  heat  power,  and  recording 
results  at  steady  state  when  the  outlet  temperature  of  the  water  did  not  fluctuate  more  than  ±  .02 
degrees  Celsius.  Successive  experiments  changed  heater  power,  then  the  procedure  was  repeated 
for  the  second  porous  sample. 

In  order  to  determine  the  heat  loss  from  the  heat  section  to  the  adjacent  section  by  conduction,  to 
the  atmosphere  by  radiation  and  convection  a  separate  experiment  was  prefonned  very  similar  to 
Pavel  [2004]. 


6.  Fluent™  Porous  Media  Results 
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Figure  3.  Fluent™  Fully  Developed  Fluid  Temperature  Profile  for  Steady-State  Porous  Media  Flow 


After  determining  the  inertia  and  viscous  terms  it  was  possible  to  set  the  boundary  conditions  in 
Fluent™.  Typical  results  for  the  fluid  temperature  within  the  porous  media  are  presented  in  Figure 
3.  Since  the  heat  was  applied  only  at  the  bottom  wall,  the  temperature  of  the  fluid  was  at  the 
highest  at  the  bottom  of  the  porous  media  and  increased  as  the  fluid  traveled  in  the  x- direction 
along  the  same  plane.  The  temperature  in  y-direction  along  the  same  plane  decreased  as  the  distance 
from  the  bottom  wall  increased.  The  Fluent™  model  was  used  to  represent  the  one-equation  (LTE) 
model. 


7.  Two-Dimensional  Numerical  Fluid  Flow  Results 

As  opposed  to  the  typical  parabola  velocity  profde  that  is  typical  of  flow  between  pipes,  the 
addition  of  the  resistance  terms  to  the  momentum  equation  that  treats  the  porous  media  as  a  pressure 
jump  resulted  in  steep  velocity  gradients  near  the  wall,  and  a  more  uniform  velocity  profile  within 
the  porous  media  region. 


8.  Laboratory  Results 

The  laboratory  experiment  mentioned  earlier  resulted  in  the  temperature  profiles  for  different 
Reynolds  numbers,  different  heat  inputs,  and  different  porosities.  The  results  below  will  compare 
the  local  heat  transfer  coefficient  with  numerical  data  based  on  the  fluid  temperature  next  to  the 
wall.  The  heat  transfer  coefficient  based  on  the  mean  fluid  temperature  will  also  be  discussed;  both 
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will  be  compared  to  the  data  obtained  from  Fluent™  data.  The  thermal  dispersion  of  the  different 
porosities  will  be  discussed.  The  porous  media  pieces  were  made  from  aluminum  with  the  thermal 
W 

conductivity  of  170 - .  The  fibre  diameter  of  the  coarse  porosity  piece  was  .004/;? ,  while  the 

mK 

diameter  of  the  finer  porosity  piece  was. 00  In? .  The  dimensions  of  the  porous  media  were  038  In?  x 
.0381/7?  x. 0127m.  The  fine  porous  media  had  a  porosity  of  .925,  while  the  coarse  porous  media 

YU  YYl 

had  a  porosity  of  .9167.  Velocities  ranging  from  .0006  — to  .002—  were  used  to  evaluate  the  heat 

s  s 

transfer  capabilities.  These  velocities  resulted  in  values  of  .6  <  Re  <  3  for  the  fine  media,  and 
.8  <  Re  <3.2  for  the  coarse  media.  The  temperature  profiles  of  the  recorded  data  showed  the 
temperatures  rising  as  the  fluid  moved  along  the  x-direction  through  the  porous  media  through  the 
interaction  with  the  heated  porous  media.  The  maximum  temperatures  were  recorded  at  the  bottom 
of  the  porous  media  due  to  the  heat  input  being  at  the  bottom.  The  temperatures  decreased  the 
further  the  fluid  was  from  the  heat  source. 


The  friction  factor  for  the  flow  in  the  fine  porosity  was  based  on  Jiang  et  al.  [1999]; 
4  PfdP  A P 


f  = 


l-£  3  M2  L 


(22) 


Where,  the  mass  flu x,M  ,  is  M  =  p fu  aiuUy  =  - 


Vol-Vol 


Vol, 


Jiang  et  al.  [1999]  also  presented  a  correlation  for  the  friction  factor  as  a  function  of  Reynolds 
number 


,  36.4 

/  = - +  .45 

Re 


(23) 


This  value  is  valid  for  Ree<2000.  It  was  found  that  the  experimental  data  correlated  well  with 
equation  23. 


The  heat  transfer  coefficient  was  found  using  the  entrance  temperature  of  the  fluid  as  recorded  by 
thermocouple.  The  results  can  be  seen  in  the  figures  below.  As  mentioned  by  other  authors,  the 
entrance  effects  were  hard  to  capture  with  the  laboratory  experiment.  The  Fluent™  line  represents 
the  average  of  the  results  of  the  different  Reynolds  numbers  for  specific  heat  input.  Typical  graphs 
are  represented  below.  The  graphs  show  that  the  heat  transfer  coefficient  was  at  a  maximum  in  the 
beginning  of  the  flow  channel.  This  was  due  to  the  mixing  that  was  occurring  at  the  entrance  of  the 
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porous  media.  The  numerical  model  captured  this  well,  while  the  laboratory  model  was  not  as 
accurate.  The  Fluent™  model  over  predicted  the  heat  transfer  coefficient. 


x/De 


—■—Fluent  ■  Re  =  .83  ■  Re  =1.47  ■  Re  =  2.37  ■  Re  =  3.17  a  Num(Re=.83)  a  Num(Re=1.47)  a  Num(Re=2.37)  a  Num(Re=3.17) 


Figure  5.  Heat  Transfer  Coefficient  Based  on  Entrance  Fluid  Temperature  at  1.64kW/mA2  Coarse  Media 


Figure  6.  Heat  Transfer  Coefficient  Based  on  Entrance  Fluid  Temperature  at  1.0kW/mA2  Fine  Media 


The  amount  of  surface  area  in  porous  media  was  large;  therefore  the  amount  of  water  in  contact 
with  the  porous  surface  was  also  large.  The  amount  of  heat  transfer  credited  to  the  thermal 
dispersion  given  in  equation  (17)  increased  as  the  Reynolds  number  increased.  The  thermal 
dispersion  for  each  of  the  Reynolds  number  is  plotted  below.  From  the  graph  it  can  be  seen  that  the 
amount  of  heat  transfer  due  to  thermal  dispersion  is  small  at  the  low  Reynolds  number.  Therefore, 
it  can  be  neglected. 
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9.  Verification  of  the  Local  Thermal  Non-Equilibrium  Model 


In  order  to  validate  the  use  of  the  LTNE  model  different  approaches  were  used  to  quantify  the 
resulting  data.  Whitaker  [1999]  proposed  methods  for  the  validation  of  the  use  of  the  two-equation 
model  based  on  the  volume  averaging  technique. 

In  order  to  validate  the  use  of  only  the  one-equation  model  Whitaker  proposed,  through  rigorous 
mathematical  modelling,  that  if  certain  conditions  were  met  the  use  of  a  one-equation  model  was 
valid. 

Whitaker’s  restrictions  to  use  of  the  thermal  equilibrium  model  were  as  follows: 

1 .  Either  s g  or  sa  tends  to  zero. 

2.  The  difference  in  the  (5  -phase  and  a  -phase  physical  properties  tends  to  zero. 


3.  The  square  of  the  ratio  of  length  scales, 


L 


tends  to  zero. 


The  porous  media,  when  subjected  to  the  above  criterion,  must  be  modelled  with  the  two  equation 
model.  Criteria  1  was  not  satisfied  because  the  porosity  of  the  two  phases  in  the  porous  media  did 
not  tend  to  zero.  Criteria  2  was  not  satisfied  because  the  difference  physical  properties  of  water  and 
aluminum  did  not  tend  to  zero.  Finally,  Criteria  3  was  not  satisfied  because  the  quotient  of  the 
length  between  the  phases  and  the  length  of  the  system  in  which  the  phases  were  evaluated  did  not 
tend  to  zero. 

Another  method  for  determining  whether  or  not  the  thermal  equilibrium  model  was  valid  was  set 
forth  by  Jeng  [2005].  Jeng  introduced  an  operator, Q,  that  was  a  function  of  the  ratio  of  the 
volumetric  heat  transfer,  hv ,  and  the  effective  thermal  conductivity  of  the  solid,  Aseff . 

This  was  derived  from  the  division  of  the  flow  path  of  the  heat  from  the  heated  surface  to  the  fluid 
stream  into  two  channels;  the  first  being  from  the  heated  wall  to  the  fluid,  and  the  second  being  the 
indirect  method  through  the  solid  matrix  to  the  fluid.  Also,  realizing  that  the  thermal  condition 
between  the  fluid  and  the  solid  matrix  will  be  closer  to  thermal  equilibrium  when  the  fin  efficiency 
was  smaller  a  relationship  between  the  volumetric  heat  transfer  and  the  thermal  conductivities  of  the 
fluid  and  solid  could  be  developed; 

r/h  H*  <  h  - ^ -  (24) 

/  2  -\-2  2  J 

A'seff  A'feff  "r 

Using  other  relationships  presented  in  the  paper,  Jeng’s  final  equation  was  as  follows: 


Q  = 


K 


X 


-H' 


vi  j  v 

In  order  to  satisfy  the  conditions  of  local  thermal  equilibrium  it  was  necessary  forQ  >  3  .  For  all 
other  conditions  the  two-equation  model  could  be  used.  When  the  equation  was  applied  to  the  data 


tanh 


K 


'/l 


seff 


(25) 
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obtained  from  the  laboratory  experiments  and  the  numerical  simulation  the  value  of  Q  was  below 
three. 

And,  finally,  another  method  to  ensure  the  use  of  the  two-equation  model  was  based  on  Kashan’s 
[2004]  definition  that  stated; 

In  essence,  the  LTNE  condition  is  declared  when  the  criterion  >9s.  -.9^  >.05  is  satisfied  at  any 
computational  node  inside  the  tube  domain  including  the  developing  region. 

The  LTNE  condition  was  checked  at  each  of  the  Reynolds  numbers  and  heat  input  experiments.  A 
typical  output  showing  the  difference  is  below.  As  can  be  seen,  the  value  of  .05  is  present  in  the 
computational  domain.  Therefore,  the  local  thennal  non-equilibrium  model  was  the  correct  model 
to  use. 

Furthermore,  by  investigating  the  figures  above  in  which  the  Fluent™  model  (LTE)  was  included  it 
can  be  seen  that  the  one-equation  model  did  not  match  the  data  to  an  extent  that  warranted  its  use  in 
the  thennal  modelling  of  porous  media.  The  numerical  data  and  the  laboratory  data  matched  up 
well,  which  further  justified  the  use  of  the  two  equation  model. 


10.  Conclusions 

A  thermal  non-equilibrium  model  was  developed  to  investigate  the  heat  transfer  in  porous  media.  It  was 
compared  to  the  data  obtained  from  a  laboratory  experiment  and  to  the  data  of  a  Fluent™  model.  The  results 
of  the  numerical  model  compared  well  with  the  laboratory  experiments,  while  the  Fluent™  model,  using  the 
thermal  equilibrium  equation,  over  predicted  the  heat  transfer.  The  criteria  set  forth  by  other  authors  were 
also  used  to  validate  the  thermal  non-equilibrium  model’s  usage  when  modelling  the  heat  transfer  in  porous 
media. 
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